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Abstract 


In this paper, the master equation for the coupled lossy waveguides is solved using the 
thermofield dynamics(TFD) formalism. This formalism allows the use of the underlying sym¬ 
metry algebras SU{2) and S'I7(1,1), associated with the Hamiltonian of the coupled lossy 
waveguides, to compute entanglement and decoherence as a function of time for various input 
states such as NOON states and thermal states. 


1 Introduction 

Recently, there has been a lot of interest in studying entanglement using coupled waveguides [H [2]. 
Specially designed photonic waveguides have provided a laboratory tool for analyzing coherent 
quantum phenomena and have a possible application in quantum computation [3]. The entangle¬ 
ment between waveguide modes is at the heart of many of these experiments and has been widely 
studied[3]. Using coupled silica waveguides Politi et.al.[5] have reported control of multiphoton 
entanglement directly on-chip and demonstrated integrated quantum metrology opening the way 
for new quantum technologies. They have been able to generated two and four photon NOON 
states on the chip and observed quantum interference, which further enhances the capabilities for 
quantum interference and quantum computing. Among various types of entangled states, NOON 
states are special with two orthogonal states in maximal super position thus enhancing their use 
in quantum information processing [6]. 

For the efficient use of these waveguides in the field of quantum information, the generated entan¬ 
glement should not decohere with time[7]- It is well known that losses have a substantial effect on 
the wave guides. Therefore, the time evolution of the entanglement is of interest in the context of 
quantum information processing using lossy waveguides. The entanglement between waveguide 
modes and how loss affects this entanglement has recently been studied by Rai et.al.jl], by using 
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a quantum Liouville equation. In this paper, we approach this problem from the viewpoint of 
thermofield dynamics [8j. This formalism has the advantage of solving the master equation exactly 
for both pure and mixed states, converting thermal averages into quantum mechanical expecta¬ 
tion values, by doubling the Hilbert space. The formalism thus makes it simpler to calculate the 
effects of noise and decoherence in the coupled two mode waveguide system. We look at the effect 
of different type of input states and show the efficiency of the states for quantum information 
theory. 

1.1 Thermofield dynamics 

Thermofield dynamics(TFD)[Sl El UHl [HI 112] is a finite temperature field theory. It has been 
applied to many branches of high energy physics and many-body systems. TFD has been used 
to solve the master equation, which helps in studying the temporal evolution of entanglement. In 
particular, it has been used to study entanglement in the presence of a Kerr medium using an 
SU(1,I) disentanglement theorem for arbitrary initial conditions by reducing it to a Schrodinger- 
like equation[T2l Hlj. Thus, all the techniques available to solve the Schrodinger equation can be 
used to solve the master equation. We derive the effect of losses giving rise to decoherence by 
solving the master equation exactly, using TFD, and then compute the decoherence and entan¬ 
glement properties of some two mode waveguide systems. 

The basic formalism of TFD is as follows: corresponding to the the creation and annihilation 
operators and a, which act on the physical space ?{, we introduce ’’thermal or tildian” opera¬ 
tors and a, which act on an augmented(fictitious) space ^ [ISl[l6l|T7| . The operators a and 
commute with a and a) and the sets of basis vectors {|n >} and {|h >} span the Hilbert spaces % 
and respectively, and {|n > x |m >} = |n, m > are the basis operators of the doubled Hilbert 
space. The Identity operator is \I >= We define a temperature 

dependent “ thermal vacuum state”, 

|0(/5) >= 0 >; with, G{9{(3)) =—i6{(3){aa — a^a^), (1) 

where, tanh0(/3) = and the corresponding thermal number distribution is h(/3) = — 

1). The statistical average of any operator is equal to the vacuum expectation value of the operator 
with respect to the thermal vacuum. In particular, the density operator p acting on a Hilbert 
space % is a. state vector l/o" >, ^ < a < 1 in the extended Hilbert space, so that averages of 
operators with respect to p are expressed as scalar products, 

< H >= Tr{Ap) =< >, (2) 

where ,the state \p°‘ >= /o“|/ > and p°‘ = p°‘ 0 I. Although, in some applications of TFD the 
a = 2 formalism is used, in this paper we work with the a = I formalism such that < A >< 
I\A\p >=< A\p > 

In TFD, the master equation is 




\p{t) >= 


-iH\p >, 


(3) 
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where, |/9 > is a vector in the extended Hilbert space % ® %* and —iH = i{H — H) + L, 
L is the Liouville term. Thus, the problem of solving master equation is reduced to solving 
a Schrodinger like equation. Symmetries associated with the Hamiltonian(such as SU{2) and 
SU{1, 1) symmetries) can be exploited to this equation. This makes the study of entanglement 
in lossy systems very tractable. 


2 Coupled waveguide system 


In optical communications, coupled waveguides are used as a transmission medium. The linear 
coupling between two wave guides is used to transfer the power one wave guide to another. To 
study the decoherence and entanglement properties of the two coupled waveguides, the model of 
Rai and Agarwal|lj is used. 

The Hamiltonian described by. 


H = huj{a^a + b^b) + hJ{a^b + b^a) 


(4) 


Where mode a corresponds to first wave guide and mode b corresponds to the second waveguide. 
The mode a and b obey bosonic commutation relations. The evanescent coupling in terms of 
distance between the two wave guide is given by J. The density operator has a time evolution 
given by 

—p = --[H,p] (5) 

In the presence of a damping term(system reservoir interaction), the evolution equations are 
governed by the Liouville equation 

§-p = -llH,p]+Cp, ( 6 ) 


where. 


Co = 








(7) 


where 7 is dissipation in the material of the waveguide. 

To study the entanglement and decoherence properties as function of time for the coupled wave 
guide system, we solve the exactly master equation Q using TFD. This allows us to study the 
response to the coupling of different input states, such as, number, NOON and thermal states, in 
the coupled waveguide system. In particular, in the absence of damping, an input vacuum state 
evolves into a two mode SU(2) coherent state. In the presence of damping, the vacuum state 
evolves into a two mode squeezed state and a thermal state into a thermal squeezed state. 


2.1 Two Coupled waveguides without damping (7 = 0) 

The time evolution of the density operator corresponding to the two coupled waveguides without 
damping consisting of fields in the modes a and b from eqs. ( 0 © is determined by, 

p = —ioj{a^ap + b'^bp — pa^a — pb^b) — iJ{a^bp + b^ap — pa^b — pb^a) ( 8 ) 

Now we apply the TFD-formalism(eq.Q), by doubling the Hilbert space and get the Schrodinger 
type wave equation 


\p >= -iH\p > 

here, the Hamiltonian H is given by 

H = uj{a) a + 6^6 — a^a — 6 ^ 6 ) + J{a)b + — a^h — 6 ^a) 

and can be decoupled into non-tildian and tildian parts: H = H — H, 

H = u}{a^a + b^b) + J{a^b + b^a) 

H = uj{d^a + b^b) + J{d^b + b^d) 


( 9 ) 

( 10 ) 


( 11 ) 

( 12 ) 


Then the solution of eq. @ is given by 

\p{t) >= exp[—iHt] (g) exp[—iHt])\p{0) > (13) 

where, |/o(0) > is an intial state 

It is clear from the above that, the two Hamiltonians H and H are independent. Hence, we can 
work with one of the Hamiltonians but since we are only interested in the physical states, we 
trace over the tilde states. 

To calculate the decoherence and entanglement properties the master equationQ,is solved us¬ 
ing the underlying symmetries associated with the Hamiltonians (eqns([TH [T 2 | )) . To see these 
symmetries explicitly we define the following operators, 

L+ = a^b, L- = b^a & L 3 = (14) 

which satisfy the SU ( 2 ) algebra, 

[L3,L±] = ±T± & [L+,L_] = 2L3 (15) 

with number operator, J\f = a^a + b^b. 

The Hamiltonian , in terms of the SU(2) generators, is 

H = ujM + J{L+ + L_). (16) 

The underlying symmetry of the Schrodinger like eqQ is 51/(2) (g) 517(2) and \p{t) > is given by 
\p{t) >= e“W^+-"*W^-|p(0) >; here, a{t) = iJt,. (17) 


Using the disentanglement formula [THl [T9], and taking the initial state |/ 9 ( 0 ) > as the vacuum 
state, the solution([l7|) is. 


\p{t) >= ^ ^“|p(0) > 


where, { = ((a(t)) = 

The density matrix in the number state basis is 


(18) 


N 

\Pi^)>= Cn,,naC^^,mJna,N -na><ma,N -rUa 

na,ma=0 


where. 


a 


ria^ria 


(l + |^|2)iV/2' 



and N = Ua + n^, N = nria + mi,. 


(19) 


A 




The entanglement properties are calculated by taking the partial transpose of p and computing 
the eigenvalues of the resulting matrix [ 20 ] : 


N 

\p{t)pT>= ^ -niaXma^N-na\ (20) 

na,ma=0 

^ria,naCma,ma ^ ~ ^ ~ ^ma,maC*na,na N — Ua >< ria, N — ma\- ( 21 ) 

The eigenvalues are 

N\ 


A 


Ariama — =t- 


{N-na)lna\ 
N\ 


y/iN - na)l{N - ma)lnalmal 

2.1.1 Entanglement properties of the system 


sin^"'“(Jt) cos^^ ^"'“(Jt); {for,na = ma) 

sin™’“ (Jt) cos^'^“™““"'“(Jt); {for,na 7 ^ rua 


( 22 ) 

(23) 


For a bipartite system, the entropy is defined as the von-Neumann entropy of the reduced density 
matrix traced with respect to one of the systems as 


N 


s = -Tr{palogpa) = - Ajlog Xj 


(24) 


i =0 


such that Pa = Tri,{p) and {A*} correspond to the set of eigenvalues of the reduced density 


operator. For the general number state, with the eigenvalues given by eqns(22, 23) the entropy is 


^ Ail 

^ = - Ej (iv-n)!n! 

iV' 

X {logl yv.-s-j-r ] + (2n) log[sin( Jt)] + {2N - 2n) log[cos( Jt)]} (25) 

(jV — n)\n\ 

We can also quantify the entanglement of the system by studying the time evolution for the 
logarithmic negativity [21], which is an easily computable measure of distillable entanglement. 
For a bipartite system described by the density matrix, p the log negativity is, 

En{T) = log 2 lip'll!, here, lip'll] = {2N{p) + 1) (26) 

where p^ is the partial transpose of p and the symbol |||| denotes the trace norm. Also N{p) is 
the absolute value of the sum of all the negative eigenvalues of the partial transpose of p. The 
log negativity is a non-negative quantity and a non-zero value of would mean that the state 
is entangled. 

Now, we consider various cases of optical input states: 

Case-1 : For two photon system as an input(i.e., N = 2): 

The entropy of entanglement of the two photon input state (both |1,1 > and |2,0 > ) is 
S = —4cos"^(Jt) log[cos(Jt)] — 2sin^(Jt) cos^(Jt) log[cos(Jt) sin(Jt)] — 4sin^(Jt) log[sin(Jt)] (27) 

and is shown in dotdashed curve of fig. 1(d). In this case, at time t = 0, we begin with a separable 
input state and thus the value of Entropy of entanglement(S) is zero. Then the value of ’S’ 
increases and attains a maximum value of 1.5 at Jt = 0.785212. Then decreases and eventually 











becomes equal to zero at Jt = 1.57061. Thus the state becomes disentangled at this point of 
time. At later times we see a periodic behavior and the system gets entangled and disentangled 
periodically. We do not see any interference effects. 

The logarithmic negativity of the various two photon states is different for various two photon 
states, unlike the entropy. 

Case-l(a.) If we take the input state as IV’m >= |lj 1 >, the output state is 

iV’ont >= Q!i| 0,2 > +a! 2 |l, 1 > +a 3 | 2,0 > . where, ai = —isin(2Jt)/-v/2 ,02 = cos(2Jt), 

0(3 = —z sin(2Jt)/\/2. Then we can write the log negativity entanglement of this system(from 


eq.(26)) as. 


En = log 2 [l + 2 (—i\/ 2 sin( 2 Jt) cos( 2 Jt) — sin^( 2 Jt)/ 2 )] 


(28) 


which is shown in thick curve of Fig.1(a). At time t = 0, we begin with a separable input state 
and thus the value of log negativity is zero. E]\r increases with time and attains a maximum value 
of 1.32875 for Jt = 0.42879, this is the maximally entangled state. Further, for Jt = 0.785212 we 
get the dips at Ej\f = l(the coincidence rate of the output modes of the beam splitter will drop 
to zero, when the identical input photons overlap perfectly in time), due to Hong-Ou-Mandel 
interference | 22 j and for Jt = 1.57061, E]\f vanishes. At later times, we see a periodic behavior, 
attributed to the inter-waveguide coupling(J). 

Case-l(b.) Now we take the input state as \'tpin >= |2,0 >, the output state is 
IV'ont >= /3i|0,2 > +/ 32 | 1 , 1 > +/ 33 | 2,0 > here, /3i = cos^(Jt ),/32 = -zsin(2Jt)/\/2, 

/Is = — sin^(Jt). The log negativity for this system is , 


En = log 2 [l + 2{—iV2 sin( Jt) cos^( Jt) + iV2 sin^( Jt) cos( Jt) — sin^( Jt) cos^(Jt))] (29) 


and is shown in dotted curve of Fig.1(a). In this case, Ej\[ increases and attains a maximum value 
of 1.32193 at Jt = 0.785212, then decreases and eventually becomes equal to zero at Jt = 1.57061. 
Thus the state becomes disentangled at this point of time. At later times we see a periodic be¬ 
havior and the system gets entangled and disentangled periodically. Unlike the earlier case for 
the |1,1 > input state, we do not see any interference effects. Clearly the entanglement dynamics 
of the states | 1,1 > and | 2,0 > are different. 


Case-l(c.) For two photon input NOON state: 

\ipin >= (| 2,0 > -|-| 0 , 2 >)/\/ 2 , the output state will be, I'lpout >= ai| 0 , 2 > -|-a 2 |l, 1 > +a 3 | 2 , 0 > . 
where, ai = cos(2Jt)/\/2, 02 = —z sin(2 Jt), 03 = cos(2Jt)/\/2. The logarithmic negativity of this 
state 

En = log 2 [l + z\/2sin(4Jt) -|- cos^(2Jt)] (30) 

is shown in thick curve of Fig. 1(c). At time t = 0, we begin with an entangled |2002 > state 
as input and thus the value of log negativity is one. We see periodic dips due to the Hong-Ou- 
Mandel interference, which is characteristic of NOON states. Ejy increases with time and attains 
a maximum value of 1.32875 for Jt = 0.356988, which is the point of maximal entanglement. 
Further, for Jt = 0.785212, Ejsf vanishes(we will see later that as we increase the number of 



1(a) Kb) 






Figure 1: (la) Shows the time evolution of log negativity for a two photon system: the thick 
curve shows the result for |1,1 > state and the dotted curve shows the result for the |2,0 > state, 
(lb) shows the time evolution of log negativity for a four photon system; the thick curve shows 
the result for |2, 2 > state, the dotted curve shows the result for |3,1 > state, and the thin curve 
shows the result for |4,0 > state. (Ic) shows the time evolution of log negativity for N photon 
NOON states (N=2,3,4,5). (Id) shows the time evolution of entropy of entanglement(S) for N 
photon system(N=2,3,4,5). 

photons in the NOON state, the entanglement does not vanish). At later times, we see a periodic 
behavior, attributed to the inter-waveguide coupling(J). 

Case-2 : For four photon system as an input(i.e.,A^ = 4): 

The entropy of entanglement for four photon system is 

S = —8cos®(Jt) log[cos(Jt)] — 8sin®(Jt) log[sin(Jt)] — 4sin^(Jt) cos®(Jt) log[4sin^(Jt) cos®(Jt)] 
— 6sin^(Jt) cos^(Jt) log[6cos^(Jt) sin"^(Jt)] — 4sin®(Jt) cos^(Jt) log[4sin®(Jt) cos^(Jt)]. 

(31) 

This is shown in dotted curve of fig.1(d). We begin with a separable input state(Entropy of 
entanglement(S) =0) after which the value of ’S’ increases and attains a maximum value of 2.03064 
at Jt = 0.785212, it then decreases and eventually becomes equal to zero at Jt = 1.57061. Thus 
the state becomes disentangled at this point of time. At later times we see a periodic behavior 
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and the system gets entangled and disentangled periodically. 

Now we consider the logarithmic negativity for each of the four photon states 

= {|2,2 >,|3,1 >,|1,3 >,|4,0 >,|0,4 >, and (|4,0 > +|0,4 >)/-v/2(fonr photon NOON 
state)}. 

Case-2(a.) For the input state as IV’in >= |2,2 >, the possible output states are, 

\4’out >= ai|0,4 > +a 2 | 3 ,1 > + 03 I 2 , 2 > + 04 II, 3 > +a 5 | 4 ,0 > . where, ai = —\/ 6 sin^(Jt) cos^(Jt), 
with,a 2 = cos(Jt)—sin(Jt) cos^(Jt)), 03 = cos^(Jt) 4 -sin^(Jt)—4sin^(Jt) cos^(Jt), 

04 = iY^ 6 )(sin^( Jt) cos( Jt) — sin( Jt) cos^( Jt)), and 05 = —\/ 6 sin^(Jt) cos^(Jt). 

0204 + 0205 + 0304 + 0305 + 0405 ] (32) 

is shown in thick curve of Fig. 1(b). The system starts at t = 0, with a separable input state and 
log negativity zero, which increases at Jt = 0.200242 to = 1.15181 and for Jt = 0.325477 dips 
at Ejsi = 1 due to Hong-Ou-Mandel interference [22] . The E^^ increases with time and attains a 
value of 1.77519 for Jt = 0.601094, this is the maximally entangled state, and for Jt = 0.785212 
again we get dips at Ej^ = 1.7277 due to Hong-Ou-Mandel interference. At later times, we see 
a periodic behaviorr , attribnted to the inter-waveguide coupling(J). Because of involvement of 
fonr photons, we can see the double the interference effect of the two photon system. 

Case-2(b.) For the input state as \'il)in >= |3,1 >, the possible output states are 

\i’out >= fei|0, 4 > - 1 - 62 13,1 > + 63 12, 2 > - 1 - 64 11,3 > -I- 65 I 4 , 0 > . where, hi = 2isin^( Jt) cos(Jt), 

62 = sin‘^( Jt) — 3 sin^( Jt) cos^( Jt), 63 = zy^ 6 )(sin^( Jt) cos( Jt) — sin( Jt) cos^( Jt)), 

64 = cos^(Jt) — 3 sin^( Jt) cos^( Jt), and 65 = —2z sin( Jt) cos^( Jt). 

The Entanglement for this system is, 

Etv = log 2 [l + 2(6162 + 6163 -|- 6164 -|- 6165 -|- 6263 -|- 6264 -|- 6265 -|- 6364 -|- 6365 -|- 6465 )] (33) 

and is shown in dotted cnrve of Fig. 1(b). We see the small difference in interference pattern with 
the |2,2 > state. The Hong-Ou-Mandel effect and the maximally entangled state occur at lower 
values. 

Case-2(c.) Now we take the input state as >= |4,0 >, theoutput state is 
IV'ont >= ci|0,4 > -I-C 2 I 3 ,1 > -|-C3|2,2 > -|-C4|1,3 > -I-C 5 I 4 ,0 > . here, ci = sin'‘(Jt), 

C 2 = 2i sin^( Jt) cos( Jt), C 3 = — \/6 sin^( Jt) cos^( Jt), C 4 = —2i sin( Jt) cos^( Jt), and C 5 = cos^( Jt). 
The Entanglement for this system 

Ejs! = log 2 [l -7 2 (ciC 2 -h C 1 C 3 -7 C 1 C 4 -t- C 1 C 5 -h C 2 C 3 -7 C 2 C 4 -h C 2 C 5 -h C 3 C 4 -7 C 3 C 5 -|- C 4 C 5 )] (34) 

is shown in thin curve of Eig.l(b). 

Case-2(d.) Eor four photon input NOON state: 

iV’m >= (|4,0 > -|-|0,4 >)/\/2 * 4!, the output state is, \'4>out >= c?i|0,4 > -|-(i 2 | 3 ,1 > -|-J 3 | 2,2 > 
-I-J 4 II, 3 > -I-J 5 I 4 ,0 > here, di = (sin^(Jt) -|- cos‘^( Jt))/\/2, 

^2 = i\/~(2){sm^{Jt) cos( Jt) — sin( Jt) cos^( Jt)), J 3 = —\/3sin^( Jt) cos^( Jt), 

di = iY^2)(sin^( Jt) cos(Jt) — sin(Jt) cos^(Jt)), and J 5 = (sin'^(Jt) -|- cos^(Jt))/\/2. 

The logarithmic negativity for this system, 

J/jv = log2[l -|- 2(did2 + J 1 J 3 + didi J 1 J 5 -|- d2d^ -\- £^ 2^4 + J 2 J 5 + d^di d^d^ -1- ^ 4 (^ 5 )] (35) 

s 



is shown in dotdashed curve of Fig. 1 (c). In this case , unlike for the two photon NOON state, 
the entanglement never goes to zero. This means that increasing the number of photons in a 
NOON state gives a more robust entanglement which is sustained at large times. To show this 
we calculate the entropy for entanglement and Logarithmic negativity for a 3 -photon and a five 
photon NOON state. The entropy of entanglement for three photon system is shown in thin curve 
of fig. 1 (d). For three photon input NOON state the entanglement is, 

En = log2[l + 2(6162 6163 6164 6263 -h 6264 -h 6364)], ( 36 ) 

where, 61 = (cos^(Jt) + ism^{Jt))/y/ 2 , 

62 = —fy^|)(sin^( Jf) cos( Jt) -|-i sin( Jt) cos^( Jf)), 63 = —zy^|)(sin^( Jt)) and is shown in dotted 
curve of Fig. 1 (c). Similarly for a five photon system as an input(i.e.,A^ = 5 ): he entropy of 
entanglement for five photon system is, shown in thick curve of fig. 1(d). 

For five photon input NOON state, the entanglement for this system is, 

En = log2 [1+2(6162+6163+6164-^6165-^6166+6263-^6264-^6265-^6266+6364-^6365-^6366+6465-^6466+6566)] 

( 37 ) 

61 = (cos^(Jf) + ^sm^{Jt))/^/ 2 , 62 = fY^|)(sin^(Jf) cos(Jf) — f sin(Jf) cos^(Jf)), 

^3 = cos^( Jf)—i sin^( Jt) cos^(Jt)),64 = — ■\/( 5 )(sin^( Jt) cos^{Jt)—ism^{Jt) cos^(Jt)), 

65 = iY^|)(sin^(Jt) cos(Jt) — f sin(Jt) cos^(Jt)), and cq = (cos^(Jt) — The Loga¬ 

rithmic entropy is shown in thin curve of Fig. 1 (c). We see that as the photon number increases the 
NOON state gets more and more robust and shows that ’’high-noon” states can be used for more 
precision measurements. These results are relevant in light of the recent experimental detection 
of entangled 5 -photon ’’NOON” states [53]. 


2.2 Entanglement properties for Inpnt Thermal States 


Now we consider the initial state |/9(0) > to be the two mode thermal state. Then, in the TFD 
formalism the time evolved state p{t) is 
|p(t) > = ^ e-*^i|p(o) > 




n 


N _n N 

\ ^ \ ^ 


—? Hi I I ' 

6 |na,nfe >< na,nb|6 


iHt 


(38) 


N 


N 


na=0 


ELa lELa / — 






n 


{na + 1)”“+1 

^ ^ nb =0 


rib,rib (- 


{fib + l)^b+l 


\na, N - Ua >< Ub, N - ribl (39) 


where, G{6) = —id{da — d^a^ + bb — Eb^), and Ua and Ub are thermal distribution functions . 


a 


r 


N 


'a,na ^l^|^|2^iV/2 


and N = na + rib- 


n„ 


The entanglement properties are calculated by taking the partial transpose of p and finding its 
eigenvalues, which are given by 

N\ 


(ha + l)2K+i)^(iV-na)!nJ 


sin^"^“(Jt) cos^'^ '^'^‘"{Jt)}; {for,na = Ub) (40) 


^[(^7+ I)n„+1 ^ ^(iv - n,,)!(N - Ufe)!nJufe! 




m 


sin’"'>+”“(Jt)cos2^-"'>-”“(Jt)}]; (n^ ^ Ub) 

(41) 


Q 













The corresponding von-Neumann entropy is: 


S = - g A, log A. = - (■«) (■«)1 

(fi i^"“ A/'l 

X i-°8l (a. ^ r)g.o.o., 1 <•"> 

Case-3(a.) Two photon system input(i.e., N = 2): 

The entropy of entanglement for the two photon thermal input state is 

1 1 ^71 ^ ^ 

‘^ = ~ (n+l )2 + ^os\Jt)] - ^g-^^sin2(Jt)cos2(Jt)log[^g^^sin2(Jt) cos2(Jt)] 

- /- \a 2 sin^iJt) IosI a- \a 2 sin'‘( Jt)] (43) 

(n+l)2 (n + 1)2 

and is shown in Figs.2(a,b,c). 

Case-3(b.) Four photon input(i.e., N = 4): 

The entropy of entanglement for four photon thermal input is 


(n - 1 - 1)2 

COS 

4n^ 

c,;t ,2 

(n - 1 - l)'^ 

Sin 

6# 

• 4 

(n - 1 - 1 )® 


4n® 

sin® 


■ sin®( Jt) log[ ^ p sin®( Jt)] 


^ = - (n+l)^ cos^Jt) cos®(Jt)] - ^Qg[ (^+i)io sin®( Jt)] 

A.'Y) ^ 4t7 ^ 

“ (- -L I'i4 sin^-^O C 0 S®(Jt) log[ siv?{Jt) cos®(Jt)] 

(^71^ ^71^ 

~ (- -L 1^6 sin^-^O cos4( Jt) log[ sin4( Jt) cos^( Jt)] 

[n + [n H- 1)’^ 

4??^ 477^ 

“ 1)8 sin®(>/0 cos^iJt) log[gryp^ sin®(Jt) cos2( Jt)] (44) 

and is shown in Figs.2(d,e,f). 

These figures show that for low values of ria and fib, the system sustains entanglement, but as the system 
gets more thermalized, it decoheres and the entropy of entanglement tends to zero. The thermal effect 
mimics a damping effect in the system. 


2.3 Two coupled waveguides with damping (7 ^ 0) 

We consider now, losses in the coupled waveguides due to system-reservoir interaction with ’ 7 ’ as the rate 
of loss due to the material of the waveguide. The time evolution of the density operator equation is 


\p >= -iH'\p > (45) 

with , 

H' =u;(a^o -I- — d^d — h^b) + J{a^b + b^a — a)b — b^a) 

— ij(a^a — aa + a^d^ -I- d^d -I- b^b — bb + b^P + b^b) (46) 

the following transformations. 


A + B 
72 ’ 

t At + Ht 

“ V 2 ' 

A + b 

.t A^ + 

(47) 

. -a + b 

■ 


r -A + b 

'= 72 ’ 

^ -.41 + fit 

72 ’ 

(48) 


give the Hamiltonian 

H" =uj{A'<A + B^B - il’i - H^H) -h J{-A^A + B"^B + iU - B^iB) 

- i-f{A^A + B^B + A^A + Hi'S -AA-BB + + B^&) (49) 


in 




Figure 2: For input two photon system: fig(a) shows the time evolution of entropy of entangle- 
ment(S) for different values of h, fig(b) shows entropy of entanglement(S) vs. Average number(n) 
for different values of t, fig(c) shows 3D-plot of entropy of entanglement in function of n and t; 
For input four photon system: fig(d.) Time evolution of entropy of entanglement(S) for different 
values of h, fig(e) Entropy of entanglement(S) vs. Average number(n) for different values of t, 
fig(f) shows 3D-plot of entropy of entanglement in function of h and t. 

























We diagonalise this Hamiltonian by applying a squeezing (Bogolubov) transformation mixing the real and 
tilde fields. 


D = ^iiA + vlA\ = ^ilA^ + viA, D = ^jliA +vlA^ k = ^llA^ + v^A 

E = H2B + E'f = nlB^ + 1 / 2 H, E = H2 B + v^B"^ & # = ^llB^ + V2B, 


where, 

/ri = coshri = 

^2 = coshr2 = 


(u! — J — 

- 17 - + (w - J)2|2 ’ 

Kw + J - iyp 

-ij - 1/7^ + (w + J)2|2 ’ 


vi = sinhri = 


V2 = sinhr2 = 


-\il? 


— i"f + -^72 + (w — J)^ 




- 17 + 1/72 + (w + J) 2 p 


( 50 ) 

( 51 ) 

( 52 ) 

( 53 ) 


and ri and r2 are the squeezing parameters and |/iip — = 1 and 1^*2 11^2 P = 1 - The final Hamiltonian 

is written as, 

Hf = S-\n)HAS{n) + S-\r2)HBS{r2) 

= nliD'^D) + nlib'^b) + nl{{E'^E) + nl{&E). ( 54 ) 


where. 


S'(ri) = exp[riKA^ — = ExplriA"^ — r\AA\\ 

S{r2) = ea;p[r2ifB+ — = Exp[r2B^ B^ — r^BB]; 


( 55 ) 

( 56 ) 


o _ *7 0— y/b+iBErA 17 o _ V 7 ^+(T+d )2 *7 o _ V 7 ^+(T+JP 27, 

ana iix — 2 2’^ ^2 — 2 2’^ — 2 2’^^^ — 2 2’ 

The generators of the SU( 1 , 1 ) algebra in terms of the modes A and B are given by 


<< 

II 

1 

+ 

II 

At A + At A + 1 

Ka 3 - —Y--- 

( 57 ) 

Kb- = BB, 

Kb+=B^B\ 

B^B + B^B + 1 

Kb 3 = 2 

( 58 ) 

and satisfy the commutation relations 

[Ka-,Ka+\ = 2 Ka 3 , [Ka 3 , Ka±] 

= iKA±; [Kb- 

,ATb+J = 2X33, [Kb 3 iKb±\= ±Kb±- 

( 59 ) 


The Casimir operators are 


KAo = iA^A-A^A), KBo = iB^B-B^B). 


Then the solution of eq. ( 451 becomes 




where. 


Kit) = 




, PA- =lt = pB-, 


PA 3 = - 2(7 + iJ)t = riBs, VA+ = -pt = PB+- 


( 60 ) 

( 61 ) 

( 62 ) 


By using the SU{ 1 , 1 ) disentanglement formula[T 5 ]. one can write eq.( 61 | as, 

|p(t) >={K{t)exp[TA+KA+](^xp[\n{T A3)KA3]exp[r a-Ka-] 

(g) ArB+KB+]exp[\n{TB3KB3)]exp[rB-KB-]}\p{0) > 


here. 


with 


r^± = 


2 pi± sinh 4 >i 


2 (j)i cosh (pi — pi3 sinh (pi 


& r iA — 


2pi 


2 pi cosh pi - Pa sinh pi 


l2 Ai3 




( 63 ) 

( 64 ) 

( 65 ) 


subscript i labels A, B. 




























We consider an initial state |/9(0) >= n Pm.n(0)|m, m, n, n >, in TFD notation. This gives us an exact 
solution of the density matrix : 


min(m',n') 00 / / . / ! \ ! I \ / r \ / r \ 


q'—O p'—O 
min{m,n) oo 

X E Et 

q—O p—0 


m+p—q \j n+p—q \j m \I n 


P 


P 


/ (m'+Ti'-2Q‘' + l) / (m + n-2q + l) 

X [r^+]p [r^3]^- ^^—^Pb-]^ 

^ Pm+p—q,m' +p' —q',n+p—q^n' +p' — q' (0) 


( 66 ) 


where C{t)is an overall phase factor due to K{t) = e 1 ■‘^t(l^Ao+KBo)+’^^Jt{NA+NB i)^ This is the exact 
solution for the density matrix of the coupled lossy system of waveguides. 


2.3.1 Calculation of Entanglement of the System for 7/0 

The entanglement properties are calculated by hrst taking trace of p{t) over the tilde space and then taking 
the partial transpose of p and calculating the eigenvalues of the resulting matrix similar to the case without 
damping. The eigenvalues are 

m 


Xr. 


{N - ma)\rn, 


•J sinh^™“ (0) cosh^^ {9)-, {for,ma = rih) and 9 = (V2j + t. (67) 


Xm^nty — i 


A^! 


i/(iV - ma)\{N - nb)\mjnbl 


sinh'"““''"''’( 6 *) cosh^^ (9); {for,ma ^ Ub) ( 68 ) 


Nl 

^{N - may.mj 


sinh^™" (9) cosh 


2N-2ma 


( 0 )] 


The entropy of entanglement of the the system is 

N N 

S = - '^X^logX^ = - ^ 

i—0 rria—O 

N' 

X {IosItAT -b-]] + log[sinh0] + {2N - 2ma) log[cosh6']} 

[N - maV-ma'- 

Thus, for two photon state as an input state, the entropy of entanglement of the system is 


(69) 


S = —4 cosh"*(d) log[cosh(0)] — 2sinh^(0) cosh^(0) log[2 cosh^(0) sinh^(0)] — 4sinh"^(0) log[sinh( 6 >)] (70) 
and is shown in thin curves of Figs.3(a,b,c,d). 


For a four photon state as an input state, the entropy of entanglement of the system. 


S = — 8cosh®(0) log[cosh(d)] — 8sinh®(0) log[sinh(0)] — 4sinh^(d) cosh®(0) log[4sinh^(0) cosh®(0)] 

— 6sinh^(0) cosh^(d) log[ 6 cosh'^(d) sinh‘^(0)] — 4sinh®(d) cosh^(0) log[4sinh®(0) cosh^(0)] (71) 

is shown in thick curves of figs3(a,b,c,d). When 7 goes to zero, we get the same entropy in figl(d) (without 
damping). As we increase the value of 7 , we can see the damping effect, for four photon system there is 
more damping than the two photon system. So, one can say that as we increase the input photons, the 
system will decohere more. This can also be seen by calculating the decoherence parameter. 

Since the state is Gaussian, we can use the the covariance matrix method to calculate the entanglement of 
the system by using Simon’s criterion [21]. The density matrix can be written as p'{t) = {r)RX {(j))p{t)R{<j))S{r), 
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Figure 3: The time evolution of entropy(S') for two photons(shown in thick curve) and four 
photons(shown in thin curve); fig(a) for 7 = 0 (without damping), fig(b) for 7 = 0 . 01 , fig(c) for 
7 = 0.03, fig(d) for 7 = 0.05, with J = 0.5. 


where S(r) is the squeezing matrix and R{(j)) is the rotation matrix mixing real and tilde fields. In our 


case, 6 * = 45° (see in eqs(47 48)) and ’r’ is squeezing parameter(see in eqs(52 53)), and p(0) is the ini¬ 
tial state of two mode system. Now we take the initial state / 9 ( 0 ) to be the two mode vacuum state, 
|p(0) >= |0, 0,0,0 >. To calculate the entanglement of the time evolved state p'{t) we go over to phase 
space description by following transformations. 




iPx), A"' = 




■ipx), 


1 




B = ^(y + ipy), 

Then, the covariance matrix is: 




B = ^{y + ipy), 


B’^ = ^{y-iPy)- ( 72 ) 


F(ri,r2) = 


( P 
0 
0 
0 
0 
0 
s 

V 0 


0 0 
q 0 

0 p* 
0 0 


0 

0 

0 

* 

q 

0 

t 

0 

0 


0 

0 

s 

0 

P 

0 

0 

0 


s 

0 

0 

0 

0 

0 

p 

0 


0 \ 
t 
0 
0 
0 
0 
0 


(73) 


q 




where, p = ^'^*cosh2ri, q = ^'^*cosh 2 r 2 , s = e '^"^*sinh2ri, t = e sinh 2 r 2 . 

Since the tildian fields are fictitious, we trace over them to get the covariance matrix for the physical modes. 
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4(b| 



Figure 4: (a) shows the time evolution of entanglement (-Eat) without damping (7 = 0), fig(b) shows 
the time evolution of entanglement with damping(for different values of 7 ) with r = 0.25, J = 0.5. 


V{ri,r2) = 


f p + q 0 -{s + t) 0 \ 

0 P* + q* 0 s + t 
-{s + t) 0 p + q 0 

V 0 s + t 0 P*+q* J 

The canonical form of covariance matrix is given by, 

a 7 \ 


V = 


7^ /3 


(74) 


(75) 


where, 


f p + q 0 \ / -(s + t) 0 

\Qp* + q* VO (5 + 0 


Then the separablility conditional] for any two mode state is 

DetaDetp + ~ — tr{aJ^JJ) > -{Deta + Detp) 

The symplectic eigenvalues are defined as. 


u^ = JI{A±Ja^-±} 


(76) 


(77) 


(78) 


where, A = Beta + Detp — 2Det'y = 2{p + q){p* + q*) + 2{s + t)^ 

and p = [DetV]~i = [{p + q)^{p* + q*)"^ + (p + q)^{s + t)^ - (p* + q*)^{s + t)^ - (s + t)‘^]~^ 
The entanglement of the system is 

En = max{0, —logiy-}. 


(79) 


For ri = r 2 = r, the entanglement for two mode vacuum states without damping(i.e., 7 = 0) is shown 
in fig4(a) and with damping(i.e., 7 ^ 0) is shown in fig4(b). We see that as damping increases the 
entanglement decreases, but, for low damping, the system seems to sustain entanglement to a large extent, 
so that it is quite robust for applications. 

In order to quantify the decoherence effects, we compute as 


Tr[p'^{t)] = Tr [TV < m,n\p‘^{t)\m,n >] 

m,n 


= Exp 


A’ft sinh (-\/27 + iJ)t 

(V^j + iJ)t cosh (-\/27 + iJ)t + (7 + iJt) sinh (-\/27 + iJ)t 


(80) 














Figure 5: (a) shows the time evolution of Decoherence for different values of 7 with J = 3 and 
fig(b) for J = 0.25 


The behaviour of decoherence is plotted figs5. We have considered two cases fig5(a), shows the variation 
of decoherence with time for strong coupling for various values of 7 and fig5(b), shows the evolution of 
decoherence with weak coupling. For strong coupling, the system decoheres in an oscillatory fashion and 
saturates to a non-zero value, while for weak coupling, one that for even short times, as the value of 
damping coefficient increases the system decoheres, to a very low value, very fast. 


2.3.2 Entanglement for two mode thermal state with damping 7/0 


Taking the initial state p(0) to be the two mode thermal vacuum state, the cqyariance matrix is given by. 


V{ri,r2) = 


c + d 0 e + f 0 

0 e + d* 0 -{e + f) 

e + f 0 c-l-d 0 

Vo -{e + f) 0 c*+d* / 


(81) 


where, 

c = -|- n 2 sinhfri), d = -I- n 2 Cosh^r 2 ), e = e~^'^*smh2ri 

and / = sinh 2 r 2 . 

Applying Simon’s criterion eq. we see that the system is entangled iff 

(ni -I- 712 )^ 6 “"^'’'*[cosh^ 2r — sinh^ 2r)^ + ^ > (ni ^^ 2 ) g- 27 *^(.Qg }^2 2 ^, _|_ 2 ^,^ ^g2) 

For ri = r 2 = r, and ni = n 2 = n , this condition is satisfied for the values of r given in the figure 5, 
in which we plot the logarithmic negativity as a function of n, for different values of r. We see that as 
the system not only gets less entangled for high values of 7 (quantified by r), but also for large n(external 
heat bath). So that in the presence of a heat bath the effect of damping increases and both have to be 
considered when generating entanglement in the lab by using coupled cavities. 


3 Conclusion 

In this paper, we have shown that the formalism of thermofield dynamics is a powerful tool for exact 
studies of coupled waveguide systems. Indeed, we have exactly solved the master equation associated 
with SU(2) and SU(1,1) symmetries for coupled lossy waveguides with and without damping. For coupled 
waveguides without damping, special attention has been given to the time evolution of the NOON states 
as inputs and we have shown that as we increase the photon number, the entanglement of the NOON 


Ifi 




























Figure 6: Shows entanglement (-Eat) vs. thermal distribution function(n) for different values of r 

states survives with time, thus making them extremely suitable for quantum information. The solution 
for damped systems was obtained by transforming the master equation to a Schrodinger type equation 
and applying the disentanglement formulae for SU(2) and SU(1,1). Our work extends that of Rai et. 
aim, as it gives the exact solution for the master equation, and, in addition shows how the entanglement 
behaves for input thermal states. Our results have also shown that the entanglement of the system can 
withstand a certain amount of damping, suggesting that it can be used for applications such as quantum 
computation, even if the waveguides are lossy. Furthermore we have shown the effect of an external heat 
bath on the system, by applying our methods to thermal input states. Our method shows the usefulness 
of thermofield dynamics in quantum entanglement problems, quite orthogonal to the approach given in 
ref. [25) . and allows us to handle damping in entanglement generation properties. We propose to apply this 
formalism to coupled light-atom systems, to shed further light on the effect of damping on the generation 
of entanglement. 

4 Acknowledgements 

M.N.K. wishes to acknowledge CSIR-UGC for a JRF fellowship. KVSSC acknowledges the Department of 
Science and technology. Govt of India, (fast track scheme (D. O. No: SR/FTP/PS-139/2012)) for financial 
support. We wish to thank Prof. C. Mukku and Prof. S. Chaturvedi for insightful comments. We also 
wish to thank V.Srinivasan for introducing us to the details of thermofield dynamics. 


References 

[1] D. N. Christodoulides, F. Lederer, and Y. Silberberg, Nature (London) 424, 817823 (2003). 

[2] S. Longhi, Phys. Rev. A 79, 023811 (2009). 

[3] S. Longhi, Laser an Photonic Rev. Vol.3, 243, (2009) 

[4] Amit Rai, Sumanta Das, G.S.Agarwal, Vol.18, No.6, Optical Express(2010) 6241 
arXiv:0907.2432v3 [quant-ph] and (2013), arXiv:0907.2432v3 [quant-ph] for this version. 

[5] A.Politi, M.J.Cryan, J.G.Rarity, S.Yu, J.L.O’Brien, Science 320, 646 (2008). 

[6] Lee A. Rozema, James D. Bateman, Dylan H. Mahler, et. al, Phy Rev.Lett. 112, 223602 
(2014), arXiv:[quant-ph] 1312.2012vl. 

[7] W.H.Zurek, Rev. Mod Phys. 75, 715-775(2003). 


17 













[8] H. Umezawa, H. Matsumoto and M. Tachiki, (North-Holland, Amsterdam, 1982) 

[9] L. Laplae, F. Mancini and H, Umezawa, Phys. Rev. CIO, 151, (1974) 

[10] Y. Takahashi and H. Umezawa Collect. Phenom. 2, 55 (1975); reprinted in Int. J. Mod. Phys. 
B 10, 1996, 1755, (1996). 

[11] I. Ojima, Ann. Phys. 137, 1 (1981) 

[12] H. Umezawa, Proceedings of the conf. on Thermoheld Dynamics, Banf, Canada 1993. 

[13] S. Chaturvedi and V. Srinivasan, J. Mod. Opt. 38, 777, (1991). 

[14] S. Chaturvedi and V. Srinivasan, Phy Rev A 43, 4054, (1991). 

[15] P. Shanta, S. Chaturvedi, V. Srinivasan and A.K. Kapoor, Int. J. Mod. Phys, B 10, 1573, 
(1996). 

[16] P. Shanta, S. Chaturvedi and V. Srinivasan, Mod. Phy Lett. A 72, 2381, (1986). 

[17] S. Chaturvedi, V. Srinivasan, G. S. Agarwal, J Phy A 32 1909, (1999). 

[18] A.Perelomov, “Generalized Coherent States and Their Applications”, Springer-Verlag (1986). 

[19] Kazuyuki Fujii, (2002). arXiv:quant-ph/0112090v2 

[20] G. S. Agarwal and A. Biswas, J. Opt. B: Quantum Semiclassical Opt. 7, 350 (2005). 

[21] G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 (2002). 

[22] C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. Rev. Lett. 59, 2044-2046 (1987). 

[23] I.Afek, O.Ambar, Y. Silberberg Science 328,No. 5980, 879 (2010) 

[24] R.Simon, Phy Rev Lett 84, 2726, (2000). 

[25] Y. Hashizume and M.Suzuki, Physica A , Vol. 392, Issue 17, 3518,(2013). 


